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CHARACTERIZATION OF REAL GAS PROPERTIES FOR SPACE SHUTTLE MAIN ENGINE 
FUEL TURBINE AND PERFORMANCE CALCULATIONS 


Gary J. Harloff 
Sverdrup Technology, Inc. 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

Real thermodynamic and transport properties of hydrogen, steam, the SSME 
mixture, and air are developed. The SSME mixture properties are needed for 
the analysis of the space shuttle main engine fuel turbine. The mixture con- 
ditions for the gases, except air, are presented graphically over a temperature 
range from 800 to 1200 K, and a pressure range from 1 to 500 atm. Air proper- 
ties are given over a temperature range of 320 to 500 K, which are within the 
bounds of the thermodynamics programs used, In order to provide mixture data 
which Is more easily checked (than H 2 /H 2 O). The real gas property varia- 
tion of the SSME mixture Is quantified. 

Polynomial expressions, needed for future computer analysis, for viscosity, 
Prandtl number, and thermal conductivity are given for the H 2 /H 2 O SSME fuel 
turbine mixture at a pressure of 305 atm over a range of temperatures from 
950 to 1140 K. These conditions are representative of the SSME turbine 
operation. 

Performance calculations are presented for the space shuttle main engine 
(SSME) fuel turbine. The calculations use the air equivalent concept. Progress 
towards obtaining the capability to evaluate the performance of the SSME fuel 
turbine, with the H 2 /H 2 O mixture, Is described. 


INTRODUCTION 

The SSME fuel turbopump operates with a mixture of H 2 and H 2 O at about 
330 atm and 1105 K. These conditions are such that actual gas dynamic measure- 
ments In the turbine are difficult to obtain. One computational tool at NASA 
Lewis for analyzing turbomachinery Is the combined computer programs of Katsanls 
and McNally (refs. 1 to 4). The program accepts constant thermodynamic and 
transport gas property data for each blade row. A current limitation of the 
analysis Is that the boundary layer program (ref. 2) requires the temperature 
dependence of the transport properties to be Input as a polynomial function. 

The program has thus far been limited to air properties. 

The goals of the present study were: 

(1) Determine the H 2 /H 2 O mixture thermodynamic and transport 
properties variation 

(2) Calculate the air equivalent performance 



(3) Curve fit the variation of y, C p , and Pr with temperature at the 
appropriate pressure 

(4) Calculate the H 2 /H 2 O mixture performance accounting for property 
varatlons from blade row to blade row 

(5) Modify the Katsanls and McNally computer programs to account for real 
gas effects. 

The first three goals are addressed herein. 


a 


sound speed 


NOMENCLATURE 


2 

Cact 


C 



2 

act 


C v 

H,h 

k 

M 


n 

Pr 

Pt 

P 

R Id/ R act 


R 

S 


T S 1 

Tt 


T 


W/Wcr 

Z 


actual sound speed squared, cm 2 /s 2 

Ideal/actual sound speed squared 
specific heat at constant pressure, J/g-K 
specific heat at constant volume, J/g-K 
enthalpy, J/g 

thermal conductivity, J/cm-s-K 
molecular weight 
mole fraction 
Prandtl number 
total pressure, atm 
pressure, atm 

Ideal/actual gas constant ratio 

total pressure ratio 

gas constant 

entropy, J/g-K 

sea level temperature 

total temperature, K 

temperature, K 

relative velocity/relative velocity at sonic condition 
compressibility 


2 



( dp/3p)j 

(3p/8T) p 

3 

Y 

5 

n 

x 

V 
p 

<Mj 


partial derivative of pressure with respect to density at constant 
temperature, atm-cc/g 

partial derivative of pressure with respect to temperature at 
constant density, atm/K 

relative flow angle, deg 

ratio of specific heats 

degrees of freedom 

efficiency 

transport property 

molecular viscosity, g/cm-s 

density, g/cc 

mixture property 


THERMODYNAMIC PROPERTIES 

Real gas properties for H 2 and H 2 O at elevated temperature and pres- 
sure are reported by Hendricks et al. (refs. 5 and 6), Goldberg and Haferd 
(ref. 7), and McCarty (ref. 8). All of the literature reviewed for hydrogen 
and steam present empirical relationships; theoretical relationships apparently 
do not exist. In contrast, a good theoretical real gas basis exists for air, 
see for example Hllsenrath and Klein (ref. 9). Perfect gas properties are 
available for hydrogen, Svehla (ref. 10) or the JANNAF tables (ref. 11). The 
JANNAF tables do cover real gas properties for H 2 O. Raznjevlc (ref. 12) 
presents a comprehensive collection of thermodynamic properties. 

The Roder (ref. 13) 1975 parahydrogen report covers temperatures from the 
triple point to 700 K for pressures up to 3000 atm. The valid range of temp- 
erature and pressure for the NASA computer programs GASP (ref. 5) and WASP 
(ref. 6) are given in table I. 'GASP' and 'WASP' stand for gas properties and 
water and steam properties respectively. The temperature of Interest, 1100 K, 
Is higher than the limit of 1000 K, but the programs are believed to be appro- 
priate for the current study. The computer programs GASP and WASP were modi- 
fied such that mixture properties could be easily calculated. 


MIXING RULE 

It was assumed that the gas components occupy the entire space at the 
mixture temperature and pressure. For each gas, the mass dependent properties 
were divided by molecular weight and, except for viscosity and conductivity, 
the mass Independent properties were mole fraction weighted, e.g., 


H m = nh = 


3 



~ ^S - Xn-jS-^ 
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n m1x = 2 n 1 = "* 

Pmlx = £Pl n1 

where the overbar represents mole basis properties. The other properties are 
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The GASP and WASP computer programs compute the partial derivatives 3p/3T| p 
and 3p/3p|y. Other partial derivatives are obtained by Maxwell's equations. 
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For viscosity, and thermal conductivity, the mixing rule of Wilke 
(ref. 14) was used. The mixing rule Is: 
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£ x 1 x 1 
1=1 1 1 


mix n 
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J-l 


ru 


where x^ Is either viscosity or thermal conductivity of gas 1. 


h 2 thermodynamic and transport properties 

Thermodynamic and transport properties for H 2 are presented In figures 1 
to 20 for a range of temperatures from 800 to 1200 K and pressures from 1 to 
500 atm. The real to perfect sound speed squared ratio varies from 1.0 to 
1.13, the perfect to real gas constant ratio varies from 1.0 to T.027, and 
compressibility varies from 1.0 to 1.13. The enthalpy and entropy are not 
zero at their respective reference temperatures e.g., 298 and 0 K due to the 
empirical nature of GASP and WASP. The shift Is not Important In the use of 
the enthalpy and entropy as only differences are of Interest. The quantities 
of Cn/R and C v /R approach their theoretical limits of (E + 2)/2 = 3.5 
and E/2 =2.5 where E = 5 Is the molecular degrees of freedom. 


H 2 0 THERMODYNAMIC AND TRANSPORT PROPERTIES 

Thermodynamic and transport properties for H 2 0 are presented In figures 21 
to 40 for a range of temperatures from 800 to 1200 K and pressures from 1 to 
500 atm. The real to perfect sound speed squared ratio varies from 0.65 to 
1.0, the perfect to real gas constant ratio varies from 1.0 to 7.4, and com- 
pressibility varies from 0.61 to 1.0. As previously mentioned the enthalpy 
and entropy are not zero at their respective reference temperatures. The 
quantities of C p /R and C v /R approach their theoretical limits of 
(E * 2)/2 = 4.5 and 1/2 = 3.5 where E = 7 Is the molecular degrees of 
freedom. 


. H 2 /H 2 0 THERMODYNAMIC AND TRANSPORT PROPERTIES 

Thermodynamic and transport properties for H 2 /H 2 0 are presented In fig- 
ures 41 to 60 for a range of temperatures from 800 to 1200 K and pressures from 
1 to 500 atm. The H 2 and H 2 0 mole functions of 0.8687 and 0.1313 were 
determined by a computational analysis of the combustion process (ref. 15). 
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The real to perfect sound speed squared ratio varies from 1.0 to 1.19, the 
perfect to real gas constant ratio varies from 1.0 to 1.39, and compressibility 
varies from 1.0 to 1.075. Again, the enthalpy and entropy are not zero at 
their respective reference temperatures. The ratios of Cn/R and C v /R 
approach their theoretical limits of U + 2)/2 = 3.5 and \/2 = 2.5 where 
5 = 5 Is the molecular degrees of freedom. 


The boundary layer code requires that y, Pr, and k be known polynoml- 
nals In temperature. Obtaining these polynomials removes the limitation of 
using air equivalent conditions. Thus, these transport properties for the 
hydrogen/steam mixture were fit using polynomial relationships for future 
turbine performance use. The polynomials were Initially assumed to be up to 
4th order In the temperature to sea level temperature ratio (T/T s ]) as the 
boundary layer computer program (ref. 2) needs this form. The pressure Is for 
305 atm and temperatures vary from 950 to 1140 K. The functional dependence 
of the transport properties on T/T s ] follows: 

molecular viscosity 


,< t s1 > 


0.2171 + 0.17677 



Prandtl number 


Pr = 0.62366 - 0.015834 



0.61 07 1 E— 04 



thermal conductivity 


2 

= 1.1879 + 0.032504 + 0.063622 

where: 


T , = 288.16 K, 

„(T sl ) . 0.3345E-03. 
k<T sl ) . 0.19E-02 - ^ 

An obvious problem with this approach Is that the transport properties are a 
strong function of pressure which Is not accounted for by this form of the 
equations . 

For eventual use In TS0NIC (ref. 1), MERIDL (ref. 3), and BLAYER (ref. 2) 
the options to account for real gas effects In the thermodynamic and transport 
properties Include: (1) using GASP and WASP and calling these subroutines as 

needed, or (2) curve fitting these data as functions of T and P and calling 
the empirical representation of the computer programs. The forms of GASP and 
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WASP computer programs are such that they can be use as (1) above with no 
additional work. 


AIR THERMODYNAMIC AND TRANSPORT PROPERTIES 

Thermodynamic and transport properties for air are presented In figures 61 
to 80 for a range of temperatures from 340 to 500 K and pressures from 1 to 
500 atm. Real to perfect sound speed squared ratio varies from 1.0 to 1.65, 
the perfect to real gas constant ratio varies from 1.0 to 1.85, and compressi- 
bility varies from 1.0 to 1.31. As previously discussed, the enthalpy and 
entropy are not zero at their respective reference temperatures e.g., 298 and 
0 K. The quantities of C p /R and C v /R approach their theoretical limits 
of I + 2/2 = 3.5 and 5/2 = 2.5 respectively, where 5 = 5 Is the number 
of molecular degrees of freedom. Real gas properties are given by Hansen 
(ref. 16) and Hansen and Helms (ref. 17). 


MIXTURE PROPERTIES 

For real gases, most of the thermodynamic and transport properties are 
functions of temperature and pressure. Using the static properties given by 
Jackson (ref. 18) see table II, the following property variation was 
calculated In this study: 
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From the above, the variation In the gas properties are: the gas constant Is 

about 7.5 percent lower than the Ideal gas constant, the actual speed of sound 
Is about 2 percent lower than the Ideal, viscosity decreases by about 
5 percent, and thermal conductivity decreases by 5 percent through the turbine. 
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AIR EQUIVALENT CONDITIONS 


Air equivalent conditions can be used as a performance reference, 

(ref. 19). Air equivalent conditions of the H 2 /H 2 O mixture were obtained 
for the inlet conditions. The Reynolds number and total temperature were held 
fixed between the gas mixture and the air equivalent case. The total pressure 
and hence mass flow are thus changed. The air equivalent calculation is neces- 
sitated by the boundary layer program requirement that polynomials of y, Pr, 
and k be supplied. (These polynomials for the H 2 /H 2 O mixture were pre- 
viously unknown.) Further, perfect gas assumptions are extensively used In 
the three computer programs of Interest (refs. 1 to 3) . 


Equivalent conditions are obtained when (see also refs. 18 and 19). 
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The air equivalent conditions are listed In table III and the definitions are 
given below for reference: 
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Table III lists the Rocketdyne mixture values, the previous study air 
equivalent values (ref. 20), and the present study mixture values and air 
equivalent values. 

The equivalent conditions 00 NOT account for real gas effects. They 
account for different (constant property) gases. Real gas affects can only be 
quantified by accounting for real property changes In the thermodynamic pro- 
cess. This has not yet been accomplished. 

The use of GASP and WASP provides, for the first time, a capability to 
determine real gas mixture thermodynamics and transport properties. The air 
equivalent performance limitation has been removed. 


SSME FUEL TURBINE PERFORMANCE CALCULATION 

The computer programs TSONIC, MERIDL, and BLAYER have been used to compute 
the air equivalent performance of the SSME turbine (see refs. 20 and 21). The 
TSONIC computer program solves for the 1 nvl sc Id flow In the blade to blade 
frame, the MERIDL computer program calculates the flow In the hub to tip mid- 
channel frame, and the BLAYER computer program computes the boundary layer on 
the blades and the endwall friction losses and the mixing losses. The solution 
Is quasl-three-dlmenslonal . The downstream whirl distribution and total pres- 
sure losses were adjusted, outside the programs by manual Intervention, to 
match the (suction and pressure) boundary layer edge velocity at the trailing 
edge. The TSONIC and MERIDL programs assume the gas to be Ideal and perfect, 
Invlscld, and Irrotatlonal . 

The SSME fuel turbine first stage stator has 41 blades, the first stage 
rotor has 63 blades, the second stage stator has 39 blades and the second stage 
rotor has 59 blades. The SSME fuel turbine mid channel blade geometry Is 
represented schematically In figure 81. The computed relative velocities at 
the hub, mean, and tip are shown In figure 82. The flow stagnates at the 
leading edge. At the trailing edge the suction and pressure side edge of 
boundary layer velocities meet at the trailing edge. There are no known 
experimental data for comparison. 

The loss analysis Is given In table IV for previous calculations (ref. 20) 
and the current study. The loss analysis methodology Is given In reference 4. 
As shown In table IV the largest losses for the stators Is the profile and 
mixing loss. The tip clearance loss Is the largest loss for the rotors. As 
shown In table V the first stage efficiency calculated here Is 0.902 versus 
0.893 from reference 20, and for the second stage the calculated efficiency Is 
0.915 versus 0.900 from reference 20. The aerodynamic performance Information 
Is presented In table V. Computed pressure, temperatures, flow angles, and 
velocity ratios are given for hub, mid, and shroud locations for the air 
equivalent case. 


SUMMARY OF RESULTS 

The study has developed thermodynamic and transport real gas data for the 
SSME mixture of H 2 and H 2 O. The GASP and WASP computer programs were modified 
and run together to determine mixture properties. The programs are general In 
that they can handle several gases over a range of pressure and temperature. 
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The current SSME capability at NASA Lewis Is limited to running equivalent 
air, quasl-three-dlmenslonal turbomachinery performance. An air equivalent 
case was calculated using constant real gas thermodynamic properties. This 
study has provided the tools to calculate the performance for the SSME mixture 
of H 2 /H 2 0. 

Future studies are recommended to Include variable real gas thermodynamic 
and transport properties In the TSONIC, MERIDL, and BLAYER computer programs. 
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TABLE I. - RANGES OF TEMPERATURE 
AND PRESSURE FOR GASP AND WASP 
COMPUTER PROGRAMS 
(ref., Hendricks 5 and 6) 
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TABLE II. - COMPUTEO BLADE ENTRANCE SSME TURBINE THERMODYNAMIC PROPERTIES, H 2 /H 2 0 MIXTURE 
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Rotor 1 
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a 1 053 
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.289x10-3 

.194xlO- 4 
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.187xlO- 4 

.563 
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a Rocketdyne values, letter to Povinelli (ref. 18). 


n H 2 = 0.8687 n H?0 = 0.1313 R = 376. S ft-lbf/lbm R M = 4.1164 

0/F «. 1.0364 


TABLE III. - SSME INLET CONDITIONS AND PROPERTIES, BALANCE AT 
FULL POWER LEVEL 
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376 
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28.2 

4.0 

24.9 

4.0 
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274.3 

177.7 

259.5 
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13 890 

14 541 

13 890 
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36 742 
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6918 
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1.34 
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1.34 
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BTU/ lbm R 

1.879 1 

.27 

1.940 

.27 
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1.0587 j 

1.07 : 
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.106 

.07297 
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53.35 

349.1 

53.35 

L_! 
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ARocketdyne, Feb. 1984 (ref. 18). 

“Same Reynolds number and total temperature as SSME. 
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TABLE IV. - SSHE FUEL TURBINE LOSS ANALYSIS 


This study 

Losses 

1st stage. 

1st stage. 

2nd stage. 

2nd stage. 


stator 

rotor 

stator 

rotor 

Profile and mixing 

0.012 

0.014 

0.013 

0.014 

Endwall friction 

.006 

.001 

.007 

.002 

Secondary flow 

.0050 

.008 

.003 

.008 

Incidence 

.001 

.003 

.001 

.004 

Tip clearance 

0 

.050 

0 

.034 

An blade row 

0.024 

0.076 

0.024 

0.062 

P tinlet/Psexit 

1.17 

1.18 

1.21 

1.15 

Previous study, reference 20 

Profile and mixing 

0.014 

0.016 

0.013 

0.017 

Endwall friction 

.007 

.002 

.007 

.002 

Secondary flow 

.005 

.009 

.004 

.009 

Incidence 

.001 

.007 

.001 

.013 

Tip clearance 

0 

.046 

0 

.034 

An blade row 

0.027 

0.080 

0.025 

0.075 

Ptinlet^Psexit 

1.19 

1.17 

1.24 

1.19 


TABLE V. - AERODYNAMIC PERFORMANCE INFORMATION, AIR EQUIVALENT 



— 

Hub 

Mid 

Shroud 

deg 6 

W, 

ft/ sec 

W/Wcr 

deg 6 

w, 

ft/sec 

W/Wcr 

deg 6 

w, 

ft/sec 

W/Wcr 

Stator 1 

In 

0 

286 

0.15 

0 

298 

0.15 

0 

310 

0.16 


Out 

-72.9 

1082 

.55 

-70.5 

991 

.50 

-68.3 

924 

.47 

Rotor 1 

In 

-52.2 

549 

.28 

-36.7 

453 

.23 

-20.7 

407 

.21 


Out 

+ 65.1 

926 

.48 

+65.9 

954 

.49 

+66.7 

1004 

.51 

Stator 2 

In 

34.6 

426 

.22 

33.7 

398 

.21 

33.2 

377 

: .20 


Out 

-70.9 

1172 

.61 

-69.5 

1063 

.55 

-68.2 

988 

! .51 

Rotor 2 

In 

-53.8 

637 

.34 

-42.1 

499 

.26 

-28.3 

414 

1 .22 


Out 

+61.9 

821 

.44 

+61.9 

867 

.46 

+63.5 

908 

: 


Pt, - 
lbf/ftr 


'’stage 


AP f /p f 


Exit 
swirl , 

ft 2 /sec 


Stator 1 
Rotor 1 
Stator 2 
Rotor 2 


474041 

471155 


471120 

367041 


367038 

364694 


364353 

284188 


1990 

1990 


1990 

1877 


1877 

1877 


1877 

1769 


0.994 

.779 

.994 

.779 


0.943 


.0197 

.0064 

.0101 


-401 

-93 

-425 

-50 


Btu = Cp aTt ■ 
lHm 1 


54.8; 


Btu 
TBm 2 


: Cp aT t = 52.39 Btu; Total = 107.2 Btu a 

lb~m 


Blade row 

Blade 

angle, 

deg 

Incidence, 

deg 

Exit design 
relative 
angles 6 

Deviation, 

deg 6 

Si 

0 

0 

-67.5 

-3.0 

«1 

-36.6 

.1 

64.5 

1.4 

$2 

21.0 , 

12.7 

-66.5 

3.0 

*2 

-35.6 

6.5 

60.5 

1.4 


a Rocketdyne value of work (with Hz/HjO mixture) « 314.6 Btu 

T5m 

b Positive angles are center-clockwise, see figure 81. 






Enthalpy, j/g 



Entropy, J/g-K 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Pressure, 

atm 



Figure 2. - Entropy versus temperature. 
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Figure 7. - Density versus temperature. 
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Figure 9. - Real sound speed squared 
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Figure li. - Partial derivative of pressure with respect to temperature at constant density 
versus temperature. 
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Figure 15. - Compressibility versus temperature. 
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Figure 16. - Enthalpy/(specific heat x temperature) versus temperature. 
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Figure 19. - Specific heat at constant pressure/gas constant versus temperature. 
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Figure 20. - Specific heat at constant volume/gas constant versus temperature. 
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Figure 31. - Partial derivative of pressure with respect to temperature at constant density 
versus temperature. 
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Figure 37. - Enthalpy/{gas constant x temperature) versus temperature. 
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Figure 39. - Specific heat at constant pressure/gas constant versus temperature. 
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Figure 43. - Specific heat at constant pressure versus temperature. 
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Figure 46. - Molecular viscosity versus temperature. 
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Figure 51. - Partial derivative of pressure with respect to temperature at constant density 
versus temperature. 



ORIGINAL PAGE fS 
OF POOR QUALITY 




Real/perfect sound speed squared 





Perfect/real gas constant 




Compressibility 


Pressure, 

atm 




Temperature, K 

Figure 56. - Enthalpy/tspecific heat at constant pressure x temperature) versus temperature 
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Figure 57. - Enthalpy/(ga$ constant x temperature) versus temperature. 
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Figure 59. - Specific heat at constant pressure/gas constant temperature. 





Figure 60. - Specific heat at constant 



constant versus temperature. 
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Figure 62. * Entropy versus temperature. 
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Figure 64. - Specific heat at constant volume versus temperature. 




Figure 65. - Partial derivative of pressure with respect to density at constant temperature 
versus temperature. 
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Figure 66. - Molecular viscosity versus temperature. 
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Figure 81. - Blade geometry for SSME fuel turbine mid channel. 
Positive angles are counterclockwise, see table V. 
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Figure 82. - Relative velocity at edge of boundary layer. 
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